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ABSTRACT 


Understanding the trajectory of a coronal mass ejection (CME), including any 
deflection from a radial path, and the orientation of its magnetic field is essential 
for space weather predictions. Kay et ah ( |2015b ) developed a model, Forecasting 
a CME’s Altered Trajectory (ForeCAT), of CME deflections and rotation due to 
magnetic forces, not including the effects of reconnection. ForeCAT is able to 


reproduce the deflection of observed CMEs (Kay et ah 2015a). The deflecting 


CMEs tend to show a rapid increase of their angular momentum close to the 
Sun, followed by little to no increase at farther distances. Here we quantify the 
distance at which the CME deflection is “determined,” which we dehne as the 
distance after which the background solar wind has negligible influence on the 
total deflection. We consider a wide range in CME masses and radial speeds and 
determine that the deflection and rotation of these CMEs can be well-described by 
assuming they propagate with constant angular momentum beyond 10 Rq. The 
assumption of constant angular momentum beyond 10 Rq yields underestimates 
of the total deflection at 1 AU of only 1% to 5% and underestimates of the 
rotation of 10%. Since the deflection from magnetic forces is determined by 10 
Rq, non-magnetic forces must be responsible for any observed interplanetary 
deflections or rotations where the CME has increasing angular momentum. 


Subject headings: Sun: coronal mass ejections (CMEs) 
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Introduction 


Knowing whether a coronal mass ejection (CME) will impact and the orientation of 
its magnetic held upon impact is critical for predicting space weather effects at Earth and 
throughout the heliosphere. The intensity of CME-driven geomagnetic storms (as measured 
by Dst) increases with the magnitude of the CME velocity and the southward magnetic 


held strength (Gopalswamy et ah 2008 Xiong et ah 2006a, 2007). Observations show that 


CMEs dehect, deviating from a purely radial trajectory, as well as rotate, changing their 
orientation. Coronal observations show that extreme dehections can happen close to the 


Sun (MacQueen et ah 1986 Byrne et ah 2010 Gui et ah 2011; Isavnin et al. 2014). Several 


of the dehections presented in Isavnin et al. 


exhibit latitudinal dehections exceeding 


30° below 8 Rq. The dehection motion can continue out into interplanetary space (Wang 


et al. 2004 Lugaz 2010; Davies et al. 2013; Isavnin et al. 2014; Wang et al. 2014). Some 


of these interplanetary dehections result from interactions with GMEs or other obstacles 


(Xiong et al. 2006b, 2009 Lugaz et al. 2012), however many of the observed cases have 


no such association. The exact relation between coronal and interplanetary dehections 
remains uncertain - the interplanetary GMEs may continue dehecting with constant 
angular momentum obtained in the corona or additional forces may cause acceleration at 
interplanetary distances. 

Early observations of GME dehections came from single coronagraph measurements 


and accordingly only latitudinal dehections were observed (Hildner 1977 MacQueen et al 


1986 Gremades & Bothmer 2004 Kilpua et al.|2009). These observations showed that 


during solar minimum GMEs tend to dehect toward the solar equator, but the trend 
becomes less obvious during other times of the solar cycle. More information on GMEs’ 
longitudinal behavior became possible with the advent of the STEREO spacecrafts allowing 
for use of multiple coronagraph viewpoints and geometric reconstruction techniques. These 
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observations showed that deflections occnr in longitnde as well as latitnde (Lin et ah 2010| 


Lngaz et ah 2010 Davies et ah 2013 Isavnin et ah|2013 Nieves-Chinchilla et ah 2013 


Isavnin et al.||2014 ). Isavnin et ah (2013) show that the latitndinal deflection tends to exceed 
the longitndinal deflection, however their CME sample only inclndes CMEs dnring solar 


minimnm. 


Wang et al. (2014) show that, even for a solar minimnm CME, interplanetary 


longitndinal deflections exceeding 30° can occnr. 

Dne to the large deflections in the low corona, where the magnetic pressnre exceeds 
the thermal pressnre, magnetic forces have become a popnlar explanation for the observed 


deflections. Kilpna et al. (2009) snggest that the open magnetic fields of coronal holes 


natnrally gnide CMEs down toward the eqnator. Shen et al. (2011) and Gni et al. ( 2011[ ) 
attribnte the deflection to the gradients in the backgronnd solar magnetic energy. Magnetic 
forces deflect a CME toward the region of lowest magnetic energy, which on global scales 
corresponds to streamer regions or, at farther distances, the Heliospheric Cnrrent Sheet 
(HCS). The direction of the deflection then depends on the location of the HCS thronghont 
the solar cycle. Dnring solar minimnm, when the HCS is relatively flat, primarily latitndinal 
deflections shonld occnr. As the inclination of the HCS increases deflections will occur in a 
wider variety of directions. 


Kay et al. (2013[) and Kay et al. (2015b) present a model. Forecasting a CME’s Altered 


Trajectory (ForeCAT), for magnetic CME deflections. Previous deflection models have 
simply compared the direction of the deflection and the magnetic gradients at distances 
greater than 2 Rq. ForeCAT simulates the deflection of CME from the integrated effect 
of the background forces beginning at the eruption of the CME. ForeCAT includes both 
magnetic tension and magnetic pressure gradients, and non-radial drag that opposes the 


deflection motion, although Kay et al. (2015b) And that the non-radial drag has minor 


effects unless abnormally large drag coefficients are used. To describe the CME’s radial 
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propagation, ForeCAT uses a three-phase propagation model, similar to that present in 


Zhang & Dere (2006), which implicitly includes the effects of drag in the radial direction. 


The magnetic tension force tends to be of comparable magnitude and point in the same 
direction as the magnetic gradients. ForeCAT reproduces both the global trends in the 


direction and magnitude of CME deflections (Kay et ah 2015b) and the deflection of a 


specihc observed case (Kay et ah 2015a) 


Kay et ah (2015b) found that the deflection forces quickly increase a CME’s angular 


momentum as it begins deflecting, typically toward the HCS. The magnetic forces decrease 
rapidly with distance causing the angular momentum to increase at a much smaller rate 


(see Figure 5 of Kay et ah (2015b)). For a strong magnetic background (e.g. a CME 
erupting from a declining phase active region) the initial increase in the angular momentum 
greatly exceeds any additional angular momentum gained beyond a few solar radii. For a 
weak magnetic background (e.g. a quiet Sun or solar minimum CME) the initial increase 
in the angular momentum is small and the slow continued increase at farther distances 
corresponds to a larger percentage of the total angular momentum gained by the CME by 
1 AU. 

ForeCAT can also simulate CME rotation due to the differential forces acting upon the 


CME. CME rotation is not as well observed as CME deflection but Vourlidas et ah (2011) 


infer an extreme rotation of nearly 80° below 30 Rq for the 2010 June 16 CME. In |Kay 


et ah (2015b) we present ForeCAT’s ability to determine rotation, but do not explore it in 


great detail. 

Here we focus on the evolution of the angular momentum of CMEs erupting from 
strong magnetic backgrounds. These CMEs tend to be the fastest and have the strongest 
magnetic held, leading to the more extreme space weather effects at Earth. We expand the 


work of Kay et al. (2015b) with simulations of 200 new CMEs dehecting and rotating during 
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propagation out to 1 AU. We ascertain at what distance CME deflections and rotations are 
“determined.” We want to And the distance at which the deflection forces have a negligible 
influence on the CME’s trajectory and orientation. This corresponds to the distance beyond 
which the CME propagates with constant angular momentum. 


2. ForeCAT 


In the most recent version of ForeCAT (Kay et ah 2015b) many of the simplifying 


assumptions were removed, most notably the restriction of the deflection to a single 
deflection plane. ForeCAT calculates the magnetic forces due to the solar background across 


the surface of a torus representing the CME flux rope (see Figures 1 and 2 of Kay et al 


(2015b). In ForeCAT, the deflection and the rotation are decoupled. The net force creates 


a deflection corresponding to a translational motion of the center of mass of the CME. This 
motion can be in any non-radial direction so the CME is free to deflect in three dimensions. 
ForeCAT’s one-dimensional rotation results from a net torque about the axis connecting 
the nose of the CME to the center of the Sun. Rotation changes the tilt of the CME, but 
does not affect its latitude or longitude. Observations show that CMEs can rotate in other 
directions so that the CME nose is not directed radially. The current ForeCAT model 
cannot account for these rotations, however they can be incorporated in the future through 
additional parameters describing the CME orientation. For a more thorough description of 


ForeCAT see Kay et al. (2015b). 


To describe the background magnetic held we use the Potential Field Source Surface 


(PFSS, Altschuler & Newkirk (1969) and Schatten et al. (1969)), a commonly used model 
that represents the magnetic held as the gradient of a magnetic potential. However, based 


on observations of Type II radio bursts, Mann et ah (2003) and Evans et al. (2008) suggest 


that the PFSS magnetic held may decay too rapidly with distance. The standard PFSS 
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model uses a source surface radius, Rss, of 2.5 Rq, beyond which the held is entirely radial 
and must fall as R~^ in the corotating frame. This source surface distance reproduces 


the global structure of the solar magnetic held (Hoeksema 1984), but may not accurately 
describe the ARs at low heights. We will also explore results using Rss = 3 Rq, which 
causes the magnetic held to fall less rapidly with radial distance. 


3. Results 

We simulate 100 CMEs sampling diherent masses and hnal propagation speeds. The 


SOHO/LASCO CME catalog (Gopalswamy et al. 2009) shows a range of CME masses 
between 10^^ g and 10^® g. Vourlidas et al. (2010) determine an average CME mass of 


1.6x10^® g from observations of 7668 CMEs. We simulate CMEs with masses between 10^^ 
g and 10^® g, ignoring the lowest mass CMEs as they tend to be the least relevant for 
space weather ehects. For each CME we assume a constant mass with distance, although 
coronagraph and interplanetary scintillation observations show that a CME mass can 


increase with distance due to accumulation of the background solar wind material (Howard 


et al. 2007 Vourlidas et al. 2010). The assumption of a constant mass makes these results 


upper limits on the total dehection and rotation. LASCO CMEs have speeds between a few 
tens of km s“^ to a few thousand km s“^ with an average value of 475 km s“^ Copalswamy 


et al. (2009). We consider CME speeds between 300 km s ^ and 1500 km s We simulate 


the CMEs all the way to 1 AU. 

All CMEs are initiated from the same active region (AR) in CR 2029 (April-May 2005, 
declining phase). The CMEs begin at an initial latitude and longitude of -15.4° and 17° 
and a tilt of 72° clockwise from the solar equator. The CMEs have an angular width of 


27.6° and a cross-sectional radius of 0.01 Rq {h in Figure 2 of Kay et al. (2015b)). 
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Fig. 1.— Deflected CME position at 2 Rq (a) and 1 AU (b) for the standard PFSS back¬ 
ground 2 Rq (c) and 1 AU (d) for a source surface radius of 3 Rq. See description in 
text. 
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We investigate how close to the Sun a CME’s trajectory is determined. Figure shows 
the position of the deflected CMEs at several distances. The background color contours 
represent the radial magnetic field at the surface of the Sun. The large red and blue region 
near the bottom left of each panel corresponds to the AR from which all the simulated 
CMEs erupt. The line contours indicate the radial magnetic field strength farther out. The 
CMEs can be deflected initially by imbalances in the AR magnetic held and on global scales 
the CMEs tend to dehect toward the HCS, the minimum in the background magnetic held 
intensity (black lines). Each circle represents the position of an individual CME. The size 
of the circle represents the CME mass with larger circles corresponding to more massive 
ones. The color of the circle represents the CME’s radial speed as indicated by the color 
bar. The line through each circle indicates the orientation of each CMEs’ toroidal axis. 

Figurea) and (b) show results using the standard PFSS model with i?ss = 2.5 Rq. 
Comparison between the position of the CMEs at 2 Rq and 1 AU, panels (a) and (b) 
respectively, shows that the fraction of the dehection beyond 2 Rq is a small component of 
the total amount. The CMEs initially deflect towards the west due to the gradients present 
in the AR. As the CMEs propagate outward the gradients from the AR weaken and the 
global gradients deflect the CMEs northward toward the HCS. While all CMEs in this case 
show deflection toward the HCS, none of them reach it. The slowest, lowest mass CMEs 
are the most susceptible to variations in the direction of the magnetic gradients. When the 
nearby global gradients differ signihcantly from the local gradients at the initial position 
they can change the direction of the deflection of the slowest, least massive CMEs. This 
behavior corresponds to a rapid increase in the angular momentum, followed by a decrease 
and finally a gradual increase at farther distances. 

As with the deflection, the majority of the rotation occurs below 2 Rq, and the rotation 
tends to increase with decreasing CME mass and velocity. The CMEs, initially aligned with 
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the AR polarity inversion line, rotate toward an orientation parallel to the HCS. 

Fignrej^c) and (d) show results with Rss = 3 i?©. The increase in the magnetic held 
strength causes an increase in the dehection and rotation at both distances. Despite the 
increase in the dehection, none of these CME reach the HCS. 


We seek to quantify the distance beyond which there would be a negligible change 
in the CME’s position at 1 AU if we do not include the dehection forces beyond this 
distance. We know the CME’s angular momentum rapidly increases close to the Sun, and 
the magnetic forces quickly become negligible, after which the CME dehects at a rate 
corresponding to constant angular momentum. Similarly, we can describe the continued 
rotation as the result of angular momentum conservation. 


Kay et al. (2015b) show that when a CME dehects with constant angular momentum 


its angular position as a function of radial distance, ^(R), can be described as 


«,K, . 1) ,1, 

where 6o is the angular position of the CME when it begins dehecting with constant angular 
momentum at a distance Rq and Vnj.,o and Vj. are the non-radial and radial CME speeds at 
Rq. The derivation of Equation [T] assumes that the radial speed remains constant. As r 
increases 9 asymptotes to a constant value. 


3.1. Deflection and Angular Momentum 

We determine how accurately Equation describes the deflection for diherent values 
of Ro For each Rq we determine the ratio, /, of the total deflection at 1 AU predicted by 
Equation!^ 0(215i?©), and the total simulated deflection at 1 AU, 0 sim , as 

f f)(215-Re) p 

^sim 
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Fig. 2.— The total deflection at 1 AU ((a) and (e)) and the ratio of the predicted to the 
total deflection for different values of Rq ((b)-(d) and (f) - (h)). The top and bottom rows 
correspond to source surface radii of 2.5 Rq and 3 Rq. 
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Figure]^ shows the total deflection at 1 AU and / for Rq eqnal to 2 Rq, 5 Rq, and 10 
Rq. The top row shows resnlts for the standard PFSS magnetic backgronnd. Fignrej^a) 
shows that assnming constant angnlar momentnm beyond 2 Rq leads to nnderpredicting 
the deflection less than 10% for all masses and velocities. This underprediction paritally 
comes from the small continued increase in the angular momentum and the fact the CME’s 
radial speed increases until 3 Rq, but we use the final propagation velocity for Figure]^ 
This underestimates the deflection between 2 and 3 Rq as the CME actually propagates 
slower at this distance. Figure [^c) shows that for Rq = 5 Rq the assumption of constant 
angular momentum yields underpredictions of 1% to 5%. Increasing Rq to 10 Rq has 
little effect. The largest errors occur for slow, high mass CMEs, which gain little angular 
momentum in the low corona but slowly gain more at farther distances. However, these 
CMEs exhibit deflections of less than 5° so the the underprediction of 5% corresponds to 
only 0.25°. 

Figures [^e)-(h) show the same as Figures [^a)-(d) but for i?ss = 3 Rq. The increase 
in the magnetic held strength close to the Sun causes a larger fraction of the angular 
momentum to be obtained below 2 Rq, causing / to increase for most CMEs at all 
distances. The low values of / for small masses occur when a CME’s angular momentum 
decreases as the strength of the global gradients begin to exceed the local gradients and the 
CME changes direction. For this source surface radius, the deflection of nearly all the CME 
can be described within 1% by assuming constant angular momentum beyond 5 Rq. 


3.2. Rotation and Angular Momentum 


The rotation can also be described by Equation [T] since the moment of inertia for 


rotation about the CME nose (see Kay et ah (2015b)) can be shown to be proportional to 


R^, assuming self-similar expansion. The non-radial velocity Vnr,o is replaced by the angular 
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velocity times the distance, cuoR-o- Figure shows the total rotation and / for several 
distances for both values of -Rgs, analogous to Figure 
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Fig. 3.— Same as Figure]^ but for the rotation. Note the difference in the contour range. 

The rotational angular momentum tends to noticeably increase out to farther distances 
than the angular momentum corresponding to deflection. Only 50%-80% of the total 
rotation is recovered by assuming constant angular momentum beyond 2 Rq. Assuming 
constant angular momentum beyond 10 Rq yields underestimates of the total rotation by 
10%. The larger source surface causes larger rotations but the behavior with distance is 
nearly the same for the two source surface heights. For the slowest, low mass CMEs an 
error of 10% may be significant as corresponds to 2.7° and 11°, for a i?ss of 2.5 Rq and 3 
Rq, respectively. For the slowest, high mass CMEs this error is negligible as it corresponds 
to less than 0.1°. 
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4. Discussion and Conclnsion 


The magnetic forces driving CME deflection and rotation decay rapidly with distance 
causing little acceleration beyond 2 Rq. The CME deflects beyond this distance at a rate 
corresponding to constant angular momentum, asymptotically approaching a constant 
displacement. The total simulated deflection at 1 AU can be predicted within 1% for most 
CMEs by assuming a CME propagates with constant angular momentum beyond 5 Rq. 
The rotation tends to evolve out to farther distances but can be predicted within 10% by 
assuming constant angular momentum beyond 10 Rq. We note that these distances are 
representative of the distance at which the solar wind transitions from a low to a high 
plasma f3, dehned as the ratio of the thermal to magnetic pressure. The solar wind can only 
efficiently transfer angular momentum to a CME through magnetic forces in a low plasma 


f3 environment, analogous to the transfer of angular momentum to the solar wind (Weber & 


Davis 1967[). Figure]^ shows the plasma f3 versus radial distance above the AR considered 


in this work. We use the Guhathakurta et ah (2006) density model and both versions of 
the PFSS magnetic held. ForeCAT does not require a coronal temperature so we assume 
a constant value of 3 MK, representative of the observed electron temperature above ARs 
(Sterling et al. 1997[ ). Figureshows that (5 exceeds unity above VJRq to 26Rq, with the 
distance being farther for larger source surface distances. 


CME dehection varies according to the relative positions of the HCS, ARs, coronal 
holes, and CME source region. The HCS is hat at solar minimum and warped at solar 
maximum. Throughout the solar cycle the relative importance of the local and global 
gradients may change as ARs become more numerous and stronger and the inclination 
of the HCS increases. Both factors may ahect the distance at which CME dehection is 
determined; this work has only considered a declining phase Carrington Rotation. 


While many authors have presented observations of interplanetary CME dehections 
















15 



Fig. 4.— Plasma (3 (ratio of thermal to magnetic pressure) versus radial distance for i?ss = 
2.5 Rq (black) or 3 Rq (red). The dashed line indicates (3=1. 
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they do not explicitly present the angular momentum at these distances, although it 
could be estimated from the published trajectories and white-light masses. If an observed 
interplanetary deflection has increasing angular momentum, some force must be actively 
accelerating the CME at interplanetary distances. Much of observed interplanetary 


deflection occurs in the longitudinal direction (Lugaz et ah 2010; Davies et ah 2013 Wang 


et ah 2014). The upcoming Solar Orbiter mission will reach as high as 34° heliographic 


latitude and as close as 0.28 AU heliocentric distance, providing an unprecedented view 
of longitudinal deflections. This perspective will allow for a more precise study of the 
evolution of CMEs angular momentum. 

ForeCAT includes the magnetic forces at all distances, including interplanetary space. 
Our results suggest that the interplanetary magnetic forces are not strong enough to 
influence the CMEs motion at interplanetary distances with high plasma (3 parameter. 
ForeCAT does include many simplihcations, notably the lack of enhancement of the solar 
wind magnetic held surrounding the CME due to the CME’s expansion and propagation. 
This effect will increase the magnetic dehection forces at interplanetary distances, however, 
ForeCAT’s current interplanetary forces are many orders of magnitude too small to produce 
noticeable interplanetary dehections. The magnetic forces at 50 Rq tend to be about 10“® 
their coronal values so the compressed magnetic held surrounding the interplanetary CME 
would need to be enhanced by a factor of over 300 times the ambient value. We suggest 
that interplanetary dehections at rates corresponding to increasing angular momentum 
must be accelerated by non-magnetic forces or result from the interaction of multiple CMEs 


Xiong et ah 

2006b 

2009 

Lugaz et al. 

2012) 


The interaction of CMEs with the HCS remains an important area of open research. 


The HCS can interfere with the propagation of interplanetary shocks (Odstrcil et ah 1996), 


and will likely also ahect CME propagation. None of the CMEs originating in the AR 
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considered in the work can reach the HCS, however Kay et ah (2015b) show cases where 


the CME cross nnderneath the cusp separating the streamer region from the base of the 
HCS. Whether the CME can penetrate the HCS farther out may depend on the distance of 
their interaction. Proper treatment of this interaction requires a more realistic description 
of the HCS. Observations suggest the HCS has a width of 10,000 km (6.7x10“® AU) at 1 
AU, with the surrounding plasma sheet being 30 times thicker (Winterhalter et ah [1994; 


Smith!2001). ForeCAT uses the Guhathakurta et ah (2006) density model, which creates 


a 56° wide density enhancement surrounding the location of the HCS, or a width of 1.48 


AU at 1 AU. This scale is comparable to the resolution obtained with MHD models (Opher 


et ah 2004), but ForeCAT does not capture the microphysics at the scale of the actual 


HCS. Additionally the impact of the CME upon the HCS will compress and distort the 
HCS magnetic held producing magnetic forces resisting the CME’s passage. ForeCAT 
cannot currently capture these effects, but with additional modihcations we will address 
the CME-HCS interaction. 


Since the dehection and rotation tend to be determined by 10 Rq it is essential to use 
accurate representations of the solar conditions in this distance range. Unfortunately this 
corresponds to the distance at which the current solar models are the most uncertain. The 
PFSS magnetic held model, a very commonly used model, assumes that the magnetic held 
is current-free and can be described as the gradient of a magnetic potential. The intense 
magnetic held helds of ARs, which can contribute signihcantly to the CME dehection, are 
certainly more complex than this simple current-free approximation. Additionally, the 
PFSS model tends to be driven by synoptic maps acquired over a full solar rotation and 
ARs can evolve on much shorter scales. These factors also apply to the global magnetic 
held conhguration, but tend to have less of an ehect. 

Our understanding of the solar magnetic held will greatly improve through the 
















observations by Solar Probe Plus, scheduled to launch in 2018 and reach the smallest 
perihelion of 8.5 Rq over six years later. One of the primary science goals of Solar Probe 
Plus is to “determine the structure and dynamics of the magnetic helds at the sources 
of solar wind.” Measuring the magnetic held at these close distances will help greatly 
constrain our magnetic held models. In the meantime, we suggest that the ForeCAT model 
can not only reproduce the observed dehection, but also constrain the unknown mass and 
drag coefficient as well as the background magnetic held. 


The authors thank the anonymous referee for the comments. 
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